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Abstract 

Fragmentation functions for heavy-light mesons, like the charmed D, D* 
mesons, are proposed. They rest on next-to-leading QCD Perturbative Frag- 
mentation Functions for heavy quarks, with the addition of a non-perturbative 
term describing phenomenologically the quark — ► meson transition. The cross 
section for production of large px D, D* mesons at the Tevatron is evaluated 
in this framework. 
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1 Introduction 



Much experimental and theoretical work has been recently devoted to the study of 
heavy flavour production in hadronic collisions. Theoretically, because the heavy 
quark mass is setting the scale in the perturbative expansion of QCD, acting as a 
cut-off for the infrared singularities, the most relevant features of this process are 
calculable within perturbation theory. Indeed the calculation in perturbative QCD 
of the differential and total cross sections to order has been performed |IJ, thus 
providing a firm basis for a detailed study of the properties of the bottom and charm 
quarks. 

The NLO one particle inclusive differential distribution will however contain terms 
of the kind a s ln(pr/m) which, in the large pt limit, will become large and will 
spoil the perturbative expansion of the cross section. This is reflected in a large 
sensitivity to the choice of the renormalization/f actor ization scales and hence in a 
large uncertainty in the theoretical prediction. In Ref. [0 this problem was tackled 
by introducing the technique of the Perturbative Fragmentation Functions (PFF) 
through which these terms were resummed to all orders and the cross section was 
shown to display a milder scale sensitivity. 

When considering heavy meson inclusive production, non perturbative effects are 
also quite important, especially for charmed mesons (being m c ~ 1.5 GeV), so they 
have to be estimated as better as possible for a reliable calculation of the production 
cross sections. They are normally introduced within the formalism of FF and indeed a 
next-to-leading order(NLO) analysis of FF into charmed mesons (in particular D, D*) 
including non-perturbative effects has been performed in Ref. |§ for e + e~ annihilation 
processes up to LEP energies. In this analysis however the charm component in the 
FF is considered only, the other components giving a small contribution to the e + e~ 
production cross section. 

On the other hand, the contribution of gluon-gluon and quark-gluon scattering 
subprocesses to the production cross section is relevant in hadronic collisions, and the 
gluonic component can no longer be neglected. 

Thus we proceed in this paper to the construction of a set of NLO fragmentation 
functions for D, D* mesons, including gluon, light and anticharm quark contribution. 
This set can therefore be used in the calculation of large pt inclusive production cross 
section for any hard collision process. 

2 Theoretical framework 

The general framework of this analysis is the following: we'll consider the frag- 
mentation into charm quark of any parton produced at large transverse momen- 
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turn pj> 3> m c , followed by the hadronization of the charm quark into the meson. 
Exploiting the difference in time scales of the two processes, short for the perturba- 
tive fragmentation of the parton into the charm and longer for the non-perturbative 
hadronization of the charm in a meson, we can factor the overall fragmentation func- 
tion of the parton i into the meson H in the following way: 

Dj*(z, n, m c ) = D\{z, fx, m c ) ® D*(z) (1) 

the <8> symbol meaning the usual convolution operation. This expression represents 
our ansatz for the fragmentation of a parton into a charmed meson. The non per- 
turbative part of the fragmentation is taken to be universal, i.e. independent of the 
parton which produced the charm quark via perturbative cascade. It is also inde- 
pendent of the scale at which the fragmentation function is taken: all the evolution 
effects are dumped into the perturbative part. 

The first part of the process can be calculated with purely perturbative tech- 
niques at an initial scale of the order of the charm mass, while we have to rely on 
phenomenological inputs to extract the hadronization non-perturbative effects at a 
fixed scale. Then using the DGLAP evolution equations at NLO accuracy we can 
evaluate fragmentation functions at the appropriate scale 0(pt), assuming that no 
scaling violation effects arise in the non perturbative part of the fragmentation func- 
tion. 

The calculation of the perturbative part of the process has been carried out in 
Ref . || . For the reader's convenience we shall briefly report the main results of this 
analysis. Using the factorization property, the charm quark production cross section 
in e + e~ collisions can be written as: 

dcx , . r — ^ P dcr>i , oc _^ k „ . > dz . . 

— {x, Q, m c ) = 2^ / Q, H)D i (z, ii, m c )— (2) 

dx £ J x dx z z 

where x is the energy fraction of the charm quark, Q is the center-of-mass energy 
and m c is the charm mass. Eq.(l) shows that the cross section is factorized into a 
short- distance term for the production of the massless parton i, and a parton FF 
D° into the charm quark, evaluated at a scale fi. When \x is taken to be of the order 
of m c , Df(x, fjL,m c ) is expressed in a perturbative expansion in powers of a s : 

Di(z,/j } m c ) = df\z) + ^d^iz^^nic) + 0(a 2 s ) (3) 

Then using the perturbative expansion of the l.h.s. of eq.(0) one obtains the explicit 
expression of df^ and coefficients. 

This has been explicitly done is Ref. obtaining the following set of NLO initial 
conditions in MS scheme for the fragmentation function of a charm quark, gluon and 
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light quarks respectively, into the charm quark: 



D c c (x,no) =6{l-x) + 



F 



2 



1 — x \ ° m\ 



2tt 

^ ( ^(x 2 + (l-x) 2 )log^ 
2tt m 



log^-21og(l-a;)-l 



(4) 



D c a (x, /io) = ^^(* 2 + (1 - x) 2 ) log «L (5) 



^, ft g(^^o) = (6) 

where /io is taken of the order of the charm quark mass, which we fix in our analysis at 
1.5 GeV. As obvious, D c is of order a s , and D^^ s is zero in the NLO approximation, 
being of the order of a s . Nonetheless, this component is generated at higher scales 
through the evolution with DGLAP equations, which involve a mixing of all parton 
components of FF. 

The PFF initial conditions (f|,|5|,|6|) , evolved up to the appropriate scale 0(pr) with 
NLO accuracy, can be used to evaluate the open charm production cross section in the 
large px region. Indeed this allows the resummation of potentially large logarithms of 
the kind a s \og(pT/fn c ), arising from quasi-collinear configurations, thus recovering a 
more reliable prediction of the spectrum at large transverse momentum than the 
fixed order 0(a s ) calculation. This is discussed in detail in Ref. in the case of b 
quarks and in ref. ]5| for charm photoproduction. 

In order to obtain the inclusive production of D, D* mesons, one has to take 
into account further the hadronization of the charm quark into the final charmed 
meson. The Perturbative Fragmentation Functions (§J5|,|6|) can be convoluted with a 
non-perturbative part, which we parametrize as 

D H np {z) = (n H )A{l-zYzP (7) 

- = ( (l-z) a z p dz 
A Jo 

where the parameters a, (3 and (njj) have to be extracted from comparison with 
experimental data at a fixed scale. Indeed a and (3 have been obtained by Colangelo 
and Nason in Ref. by fitting ARGUS data for D, D* mesons fragmentation 
functions in e + e~ collisions at center-of-mass energy of 10.6 GeV. They are reported in 
Table [TJ. We'll use their determination, on the ground of our universality assumption 
for the non-perturbative part of the fragmentation functions. It is worth mentioning 
that we have challenged this universality assumption by comparing our fragmentation 
functions ([[]), evolved up to 90 GeV and using the ARGUS parameters of Table [I], 
with LEP data by OPAL [0. Reasonable agreement has been found, giving support 
to our hypothesis. It is also worth noting that Colangelo and Nason original work 
gives a sub-set of our fragmentation functions only, since they were addressing 
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Meson 


a 


P 




(nil) 


D° 


1.0 


3.67 


0.58 


0.58 


D*° 


0.6 


5.4 


0.3 


0.29 




Colangelo /Nason 


HERWIG 


JETSET 



Table 1: Collection of parameters which describe the non-perturbative part of the frag- 
mentation functions. In the Colangelo-Nason paper, ref. ||, these a and (3 where obtained 
with A 5 = 100 MeV. 



the non-singlet component, by far the dominant one in the large z region. We deal 
instead with the whole set, including mixings with gluons and antiquarks: a more 
complete analysis of D mesons production in e + e~ collisions at LEP with the full set, 
extending also down to small- z values, will be presented elsewhere |§. 

The parameter (uh) in eq. (|7|) is the mean multiplicity of charmed mesons pro- 
duced in the process, i.e. how many mesons does the non-perturbative fragmentation 
of a c quark produce. In the analysis of Ref. || the normalization condition fixed by 
(uh) has not been determined. We have therefore used the mean multiplicity simu- 
lated by the HERWIG generator in the process e + e~ — > cc at 90 GeV. This provides 
us with the last missing parameter for fully defining the non perturbative part of 
the fragmentation function. The values we used are summarized in Table [l]: they 
describe either the non perturbative fragmentation of a charm into a D° or a D*°, 
or alternatively of an anticharm into a D or a D*°. As a double check, the same 
multiplicities have also been extracted from JETSET, by fragmenting a charm quark 
of 5.3 GeV energy. They are also reported in Table [I], and can be seen to agree well. 
It must be noted that the D mesons multiplicities include the feed-down from the D* 
decays. It is also worth noticing that the non-perturbative fragmentation function 
(0), with the parameters shown in Table [I], compares nicely with a Peterson fragmen- 
tation function with e = 0.06, used for instance in [E^] to describe the transition from 
c quarks to D mesons. 

The use of an explicit parametrization for the FF at a given scale extracted from 
HERWIG has been successfully made in the past for the case of the inclusive produc- 



tion of light and strange mesons [|T^]. In these papers not only the multiplicity but 
also parameters a and f3 dictating the shape of the FF (see Table were extracted 
from HERWIG at some large scale. We have refrained from doing so in this work, 
since the large mass of the charm quark provides us with a solid ground for evaluating 
in pQCD at least the perturbative part of the fragmentation function, see eq. (|I|). For 
the sake of completeness we do however provide in Fig. [I] a comparison between the 
fragmentation function which we obtain by evolving the ansatz of eq. (|l|) and that 
we get using the HERWIG parametrization given in Table 0, at a scale of 90 GeV. 
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Figure 1: Comparison between the Dj? fragmentation function pro- 
duced by HERWIG at 90 GeV and the one given by our ansatz ([[]) with 
the non-perturbative parameters by Colangelo and Nason (see table [l]) . 



The value of A 5 = 100 MeV has been used in the evolution, for consistency with the 
Colangelo-Nason fits jj. The Sudakov form factor j| has also been included, as it 
was taken into account in ref. Its effects are however small, of the order of a few 
percent, both on the fragmentation function itself and on the cross sections which 
will follow. 

A number of comments about this plot are in order. The discrepancy in the small 
x region is due to the PFF being enhanced by mixing with the gluon splitting kernel, 
while the HERWIG result is suppressed by phase space constraints. This difference is 
however of no practical importance in the evaluation of hadronic cross sections, since 
such small x values probe the very large px tail of the kernel cross section, where the 
latter is small. The discrepancy around the maximum can instead be considered as 
a normalization off-set, which could be eliminated by a fine tuning of both FFs to 
some experimental data. It is however worth noticing that the present accuracy of 
the data is not better than the uncertainty originating from the difference of the two 
FFs. Moreover, the consequences on the observable hadronic cross section are small. 
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a 










2.75 ± 0.08 
1.12 ± 0.51 


2.72 ± 0.03 
0.36 ± 0.19 


53.6 
0.0045 


0.58 
0.0013 


c^D*° 


2.01 ± 0.06 
0.16 ± 0.08 


2.42 ± 0.08 
0.016 ± 0.03 


12.4 
0.008 


0.30 
0.007 



Table 2: Collection of parameters which describe the fragmentation functions produced 
by HERWIG at 90 GeV, parametrized as D?(z) = N l H (l - z) a 'z^\ It also holds (n) = 
NB(a + 1, (3 + 1), B being the Euler Beta function representing the normalization integral 
of the FF. Note that the parameters a and (3 do not have to coincide with those of eq. @), 
because they include also the perturbative component of the FF. 



3 Results 

Using the perturbative initial conditions (§J^J^) and the non-perturbative parametriza- 
tion (0) with the parameters summarized in Table [I], we have a NLO evaluation of 
the parton FF in D, D* mesons at a scale fio ~ m c . Then using the DGLAP evo- 
lution equations to NLO accuracy one can evaluate the fragmentation functions set 
at any desired factorization scale fi. By convoluting this set with the NLO kernel 
cross sections for massless parton scattering |TT[] one gets a prediction for D and D* 



production at large px at a hadron collider. Figures ||a,b show the px and rapidity 
spectra for the Tevatron. Also shown (solid line) is the pure charm quark NLO cross 
section, as predicted by the Perturbative Fragmentation Function approach. These 



results have been obtained with the MRS-A structure functions set |14| and with a 
A 5 value of 100 MeV. The Sudakov form factor || has also been included: we notice 
here once more that its effects are numerically small, below 10%, in the pr-rapidity 
range we have considered. 

As expected the D mesons cross sections lie below the charm quark one. It can 
however be checked that summing them up and introducing an additional factor of 
two to allow for D + and D* + state (which are assumed to fragment like the neutral 
ones) the charm cross section is almost reproduced (wide-dotted line). The small 
residual gap is given by the mesons FF being softer and therefore producing a smaller 
cross section. 

Many uncertainties do of course affect this result, though not shown in the plots. 
First of all, the mesons fragmentation functions will share all the uncertainties related 
to the heavy quark PFFs on which they are built on. Factorization scale and initial 
scale dependencies, of the kind studied in 0, will also appear here with a similar 
behaviour, leading to an uncertainty of order 20-30%. These fragmentation functions 
will also share the same shortcomings of the PFFs. This means that their description 



6 





D° 


! \\ 


D*° 




charm 




mesons sum 


I MRS-A 




r Vs= 1800 GeV 




: y = 

i t i t 1 i i i t 1 i i t t T i t t t T i t t i 1 





MRS-A D 

Vs= 1800 GeV D 

p T = 20 GeV charm 



P T (GeV) 

Figure 2: Cross sections for mesons production at the Tevatron, as 
predicted by our fragmentation functions approach. Comparison with 
the pure charm quark cross section is also shown. The mesons sum also 
includes D + and D* + contributions. 



is not accurate at low pt and at the edges of phase space (see also || for a discussion of 
this point), where unresummed higher order corrections and non-perturbative effects 
play an important role. This leads to the impossibility of performing a meaningful 
comparison with the D* photoproduction data collected at HERA, since the minimum 
Pt, of the order of 2-3 GeV, is too low and the edges of phase space in rapidity are 
probed. 

A second kind of uncertainty is related to the determination of the non-perturbative 
parameters a, (3 and (nn)- The set we have chosen was fitted to ARGUS data, and 
has been picked mainly for illustrative purposes, although, as stated above, it re- 
produces fairly well the actual OPAL data from LEP. An analysis of all LEP data, 
when available, will certainly lead to a more precise determination of these parame- 
ters. Of course also a detailed measurement of the D, D* inclusive cross sections at 
the Tevatron would be very helpful and give complementary information on the FF, 
particularly for the gluon terms. 

To conclude, we have presented a model for the D and D* fragmentation func- 
tions based on PFF for heavy quarks complemented by a factorized non-perturbative 
term describing the quark-meson transition. Predictions have been given for large pr 
charmed mesons production at the Tevatron. 
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